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INTRODUCTION 

Controlled  rectifiers  employing  power  semiconductors  are  used  to  vary  the  average 
power  supplied  by  a  constant  voltage  a.c.  source  to  a  load  circuit.  One  of  the  problems 
inherent  in  the  design  of  such  a  system  is  the  determination  of  the  commutation  or 
overlap  angle.  Knowledge  of  the  commutation  angle  gives  the  designer  information 
about  the  rate  of  change  of  current  in  the  power  semiconductor  and  allows  him  to 
accurately  model  the  load  current. 

A  short  review  of  the  literature  failed  to  produce  a  rigorous  determination  of  the 
commutation  angle.  Common  methods  of  analysis  either  ignore  the  commutation 
reactance  contained  in  any  rectifier  circuit  or  assume  constant  load  current. 

The  model  developed  here  represents  the  circuits  depicted  in  figures  1  and  2.  A  full 
wave  rectifier  employing  a  transformer  connected  across  the  single-phase  source  and  a 
bridge-rectifier  circuit  are  shown  in  figure  1 .  The  circuit  shown  in  figure  2  also  gives  full 
wave  rectification  and  consists  of  a  center-tapped  transformer  and  controlled  rectifiers. 
While  the  center-tapped  transformer  has  the  effect  of  transforming  the  single-phase 
source  into  a  two-phase  source,  the  circuit  is  commonly  referred  to  as  a  single-phase 

transformer.1 

The  load  current  obtained  from  either  circuit  is  the  same;  the  choice  of  which  circuit 
should  be  employed  is  application  dependant.  Factors  to  be  considered  in  such  a 

selection  include  cost,  required  load  voltage,  and  the  available  voltage  supply.1 


DESCRIPTION  OF  MODEL  AND  BOUNDARY  CONDITIONS 

The  method  used  to  calculate  the  commutation  angle  and  load  current  is  dictated  by 
the  fact  that  no  a  priori  knowledge  of  the  magnitude  of  the  current  at  the  boundaries 
exits.  The  current  at  the  boundary  points  depends  on  the  firing  and  commutation 
angles  as  well  as  the  circuit  parameters.  As  a  result,  the  conventional  Laplace  transfor¬ 
mation  approach  is  not  applicable. 


’Dewan,  S.  B.  and  Straughen,  A.,  Power  Semiconductor  Circuits,  John  Wiley  and  Sons. 
New  York,  1975,  p  216.  ~ 


To  solve  the  problem,  a  set  of  equations  with  undetermined  constants  of  integration 
were  developed.  The  equations  represent  the  current  in  the  load  during  commutation 
and  conduction  and  the  phase  "a"  current  during  commutation.  Through  successive 
application  of  the  boundary  conditions,  four  independent  transcendental  equations 
containing  the  constants  of  integration  and  commutation  angle  were  generated.  The 
equations  were  then  solved  simultaneously  through  an  iterative  process  resulting  in 
numerical  values  valid  for  the  input  circuit  values  and  firing  angle. 

Continuity  of  load  current  was  assumed.  The  method  of  application  of  the  boundary 
conditions  forces  continuous  load  current  and  cases  resulting  in  a  discontinuous  re¬ 
sponse  will  not  lead  to  accurate  results.  Fortunately,  it  is  not  difficult  to  determine  if  a 
discontinuous  current  situation  exists.  The  procedure  to  determine  continuity  of  current 
is  discussed  later. 

During  the  conduction  phase  of  operation,  either  circuit  may  be  represented  by  the 
single  phase  equivalent  shown  in  figure  3.  Since  the  response  of  both  "phases”  to  the 
load  is  symmetrical,  it  is  only  required  to  calculate  the  phase  "a"  response  for  one  half 
cycle.  The  response  during  the  second  half  cycle  is  simply  a  time  shifted  duplication  of 
the  phase  "a”  current. 

The  commutation  circuit  (fig.  4)  represents  the  circuit  while  SCR  1  is  coming  on  and 
SCR  2  is  going  off.  The  voltage  sources  are  sinusoidal  so  vg  =  Vsincot  and  vb  =  Vsintot. 

The  time  reference  is  established  so  that  t  =  0  when  the  phase  "a"  voltage  becomes 
positive.  SCR  1  is  fired  at  cot  =  a.  It  is  assumed  that  SCR  1  is  forward  biased  at  this 
point.  During  commutation,  both  phases  conduct  simultaneously  and  the  load  current. 

id,  is  the  sum  of  the  phase  currents  ig  and  ifa. 

The  expected  load  response  is  illustrated  in  figure  5.  The  following  boundary 
conditions  used  in  the  formulation  of  the  load  current  equations  are  readily  apparent 
from  reference  to  the  figure: 

ia(«)  =  0 

id(o.  +  p)  =  ia  (a  +  p) 
id  (a  +  p)  =  iac  (a  +  p) 
iac  (rx  +  n)  -  id  (rx) 


■  *■  -r»'^  f.*  nj»  .TXJUl  JunuruxK*  it*  vi  vn  TWVSW^JWDWWIWl  v.  v  x.-  \  »„  wit  »ytv  \r;  w  v_-«-.  v  -,*  -  <-.- 


DEVELOPMENT  OF  EQUATIONS  FOR  LOAD  CURRENT  AND  CONSTANTS 


Beginning  with  the  method  described  by  Hoft  ,  the  voltage  balance  around  the 
outside  and  lower  paths  in  figure  4  is  written 


-Vs in  ox  +  L,  -jj-  +  L 


d  ia  i  ti(ia  +  ifj ) 


+  R(ia  +  i),)  +  Vc  =  0 


Vsinttit  +  Lt  +  L  -  a^-~—  +  R (  C  +  ‘b )  +  I'r  =  0 

Adding  these  two  equations  and  defining  a  load  current  during  commutation, 
id  =  ia  +  ib,  leads  to  an  equation  without  a  forcing  term 

L  +  _^=o 

clt  +  L,  +  2L  +  L,  +  2L 

Equation  1  is  a  linear  first  order  differential  equation  and  has  the  solution 


.  -  l.  -  <*)/0)T,  _  _< 

‘(I  K  I  1  D 


where 


L,  +  2  L 


and  k,  is  to  be  determined  from  the  boundary  conditions.  t1  is  the  time  constant  of  the 
circuit  during  commutation. 

In  order  to  obtain  an  independent  expression  for  the  phase  "a"  current  during 
commutation,  consider  the  voltage  balance  around  the  outer  loop  of  figure  4 


-VsitHMt  +  —7-  +  \  n  =0 

at 


where 


'  i)  -  L  (j‘  +  R  i,i  +  ' 


^Hoft,  R.G  .  Semiconductor  Power  Electronics,  Van  Nostrand  Reinhold  Co.,  New  York, 
1986,  p  126. 
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Taking  the  derivative  of  equation  2  and  substituting  for  VQ  and  id  in  equation  3  yields  a 
separable  differential  equation  for  i  which  can  be  integrated  as 


'•te 


j-  sintat  -  jjr  e~Um~a),mx'  idt 


-ir  cost* +  4  e  +k2 

CO  2. 


To  determine  the  constant  k2,  apply  the  first  boundary  condition  (ig(a)  =  0)  to  equa¬ 
tion  4.  Then, 

(5) 

Applying  the  second  boundary  condition  [id(a  +  p)  =  ig(a  +  p)]  in  equations  2  and  4  and 
substituting  for  k2  from  above  results  in  an  expression  for  k1 


2  f_L 

1  +  i  o)L, 


k\  =  T~T~u 7^7  ™ r  l(wa  -  +  -5- 


Now  consider  the  conduction  period  for  phase  "a”  (a  +  p  <  cot  <  a  +  n).  Referring  to 
figure  3,  the  voltage  around  the  circuit  can  be  expressed  as 

—l  sin  on  +  (Li  +  L)  ■ -  +  R  \ ,  =  0 

The  solution  for  the  conduction  current  is 


*  .  .  „  _K.  ,  ,  How  -  -  iD/on-* 

=  y  Sin  (Ml  -  <P)  +  ku’  •  -  R 


where 


7.  =  V ( to/- ):  +  R~ .  <t>  =  tan 


I,  +  /. 


,  <')</.  +  /,, ) 


t2  is  referred  to  as  the  conduction  time  constant  and  k3  is  to  be  determined  from  the 
third  boundary  condition. 

The  third  boundary  condition  deserves  some  explanation.  Note  that  id(a  +  p)  = 
ia(a  +  p)  +  ib(a  +  p)  and  ib(a  +  p)  =  0.  In  addition,  ig(a  +  p)  =  iac(a  +  p);  therefore, 
id(a  +  p)  =  igc(a  +  p).  Substituting  this  condition  into  equations  7  and  2  leads  to  an 
expression  for  k3 

k\  =k\  c~M/WTl  -  ^  sin  (a  +  p  -  d>)  (8) 

One  more  independent  equation  is  required  to  completely  specify  the  four  un¬ 
knowns  (kr  k2,  k3,  and  p).  To  obtain  the  last  equation,  the  remaining  boundary  condi¬ 
tion  is  used.  Note  first  that  for  steady  state  operation,  igc(«  +  n)  =  i b(u).  Using  reason¬ 
ing  similar  to  that  above,  since  ia(ot)  =  0,  igc(a  +  n)  =  id(a).  Using  this  information  in 
equations  2  and  4  yields 

k\=-j  sin  (a  -  O)  +  ky  ir(K  “  ^,/wt2  (9) 

Equations  5,  6,  8,  and  9  are  four  independent  equations  defining  kr  k2,  k3,  and  p. 

Substituting  equations  6  and  8  into  equation  9  yields  a  single  (albeit  cumbersome) 
transcendental  equation  for  p  in  terms  of  circuit  parameters.  Numerical  solution  to  this 
equation  then  allows  calculation  of  kv  k2,  and  k3.  Once  these  constants  are  known,  igc, 

id,  and  ig  can  be  readily  calculated. 

DETERMINATION  OF  MINIMUM  PHASE  CONTROL  ANGLE 

As  mentioned  previously,  it  is  assumed  that  SCR  1  is  forward  biased  at  tot  =  a. 
Once  p  is  numerically  determined,  it  is  a  simple  matter  to  check  this  assumption.  Just 
prior  to  the  instant  SCR  1  is  fired,  Kirchoff's  voltage  law  can  be  written  for  the  top  loop  in 
figure  4 


Since  ia  =  0,  ib  =  id  and,  using  equation  2,  the  voltage  on  the  SCR  is 


L,k\ 

Vak\  =~JrL  +2Vsin<at 


If  VAK1  is  less  than  zero,  then  the  assumption  was  violated  and  the  SCR  could  not  be 


forwarded  biased  at  cot  =  a.' 


DETERMINATION  OF  CONTINUITY  OF  CURRENT 


In  the  previous  derivations,  it  was  presumed  that  operation  was  under  steady-state 
continuous-current  conditions  (i.e. ,  id  is  never  zero).  There  are  two  indicators  which  can 

be  examined  to  determine  if  this  assumption  is  violated. 


Some  combinations  of  circuit  parameters  and  firing  angle  always  produce  a  discon¬ 
tinuous  current.  These  can  be  determined  by  assuming  a  discontinuous  current  (so 
there  is  no  commutation)  and  writing  Kirchoff's  law  for  phase  "a"  conduction 


Vsinon  =  (Lf  +  L)  +  Ria  +  Vc 


This  equation  applies  for  a<cot<P  where  p  is  the  electrical  angle  at  which  i  returns  to 
zero.  Solving  the  voltage  equation  for  i  and  writing  i  (cot  =  P)  produces  a  transcenden¬ 


tal  equation  for  p. 


£-  jsin  (a  -  R(a  P)/wL  +  ^  .v/«(p  -  <J>)  -  ^  =0 


Hoft,  R.G.,  Semiconductor  Power  Electronics,  Van  Nostrand  Reinhold  Co.,  New  York, 
1986,  pi  26. 


Dewan,  S.B.  and  Straughen,  A.,  Power  Semiconductor  Circuits,  John  Wiley  and  Sons, 
New  York,  1975,  chap.  3. 
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Since  ia  is  always  positive  while  SCR  1  conducts,  it  is  not  necessary  to  solve  this 

equation.  Subsitution  of  values  for  (J  which  are  less  than  the  true  value  will  produce  a 
positive  number  from  the  above  equation,  while  substitution  of  values  for  (3  which  are 
greater  than  the  true  value  will  produce  a  negative  number.  If  substitution  of  the  value 
(3  =  a  +  n  results  in  a  negative  value,  then  the  true  (3  is  less  than  a  +  n  and  the  circuit 
shuts  off  completely  prior  to  firing  SCR  2.  In  this  case,  the  current  is  discontinuous. 

The  above  criterion  for  discontinuity  is  not  sufficient  in  all  cases.  Some  cases  which 
pass  that  test  still  result  in  negative  values  for  phase  "a"  current  (indicating  a  discon¬ 
tinuous  case)  for  some  value  of  cot  between  a  and  a  +  n.  These  cases  can  only  be 
found  by  calculating  ji  and  the  circuit  response  for  the  range  of  values  a<cot<a  +  rr. 


RESULTS 

The  transcendental  equation  for  p  was  solved  numerically  on  a  CDC  Dual  Cyber 
170/750.  The  program  (Appendices  A  and  B)  allowed  input  of  all  circuit  parameters  and 
a  f:ring  angle,  a.  The  commutation  angle,  p,  was  then  calculated  using  a  straightfor¬ 
ward  method  of  bisections.  If  the  input  a  failed  the  check,  as  discussed  previously,  a 
was  incremented  by  1  degree  and  p  was  recalculated.  This  continued  until  an  accept¬ 
able  a  was  found  (to  within  1  degree).  The  current  values  were  then  calculated  and 
plotted.  In  spite  of  the  unsophisticated  method  of  bisections,  run  times  were  relatively 
short. 

The  plots  from  several  runs  are  included  which  show  some  trends  as  various 
parameters  are  varied  (figs.  6  through  11).  The  input  a  for  all  runs  was  zero  so  that 
the  «  indicated  was  the  minimum  a  determined  iteratively.  In  general,  large  values  of 
inductances  were  selected  to  clearly  display  various  trends. 

Most  runs  provided  current  plots  with  the  expected  shapes,  even  for  the  flat  re¬ 
sponse  expected  from  a  very  large  load  inductance  (fig.  10).  One  run  in  which  the  load 
electromotive  force  is  50  volts  (fig.  1 1 )  shows  a  negative  phase  "a"  current  indicating  the 
discontinuous  case,  which  violate  our  initial  assumptions.  Therefore,  these  results  are 
invalid.  This  case  is  interesting  because  it  did  not  have  a  value  of  (3  greater  than  a  +  n. 
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P K  SOURCE  VOLTACE  (V) 
SOURCE  f REOUENCY  (HZ) 
fIRINC  ANCLE  (DEC) 
TRANS  INDUCTANCE  (H) 
LOAD  INDUCTANCE  (H) 


100  00  LOAD  RESISTANCE  (OHU)  100  00 

eo  oo  load  Etr  (v)  so  oo 

10  00  CONDUCTION  TAU  (S)  0  0060 

0.3000000  COLMJTATION  TAU  (S)  0.00*5 

0.3000000  COUAJTATION  ANG  (DEC)  2S.3335 


Figure  9. 


Load  response  for 


back  EMF  increased  to  30  V 


PK  SOURCE  VU.TACE  (V) 
SOURCE  FREQUENCY  (HZ) 
riRINC  ANCLE  (DEC) 
TRANS  INDUCTANCE  (H) 
LOAO  INDUCTANCE  (H) 


100. 00 
60  00 
1. 00 

0  2500000 
100. 000000 


LOAD  RESISTANCE  ( OHM)  100  00 

LOAD  EUT  (V)  0  00 

CONDUCTION  TAU  (S)  1  0025 

CONFUTATION  TAU  (S)  1  OOU 

CONFUTATION  ANC  (DEC)  56  0005 


Figure  10  . 


Flat  response  for  large  load 


inductance 


I  3 


_  phase  a  current 

_ LOAD  C04A4UT AT  ION  CURRENT 


•J.M  l«l  00 

wr  (oec) 


PK  SOURCE  VOLTAGE  (V) 

100. 00 

LOAD  RESISTANCE  (OKU) 

100. 00 

SOURCE  FREQUENCY  (HZ) 

SO.  00 

LOAO  EIT  (V) 

50.00 

FIRING  ANCLE  (OEG) 

11. 00 

CONDUCTION  tau  (S) 

0. 0060 

TRANS  INOUCTANCE  (H) 

0.  3000000 

COMMUTATION  TAU  (S) 

0.  0045 

LOAO  INOUCTANCE  (H) 

0.  3000000 

COMMUTATION  ANC  (DEC) 

S. 7505 

Figure  11.  Invalid  response  showing  negative  current  swing 

(Current  in  this  case  is  discontinuous  and  is  not 
correctly  calculated  with  this  model.) 
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FORTRAN  SOURCE  CODE 
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PROGRAM  SGLFAZ ( INPUT, OUTPUT, TAPE5=INPUT, TAPE6=OUTPUT, TAPE7 ) 

A***************************************************** 

* 

*  FORTRAN  PROGRAM  SINGLE  PHASE 


*  JOHN  PAPPAS  &  JOSEPH  BENO 

*  UNIVERSITY  OF  TEXAS  AT  AUSTIN 

*  DEPARTMENT  OF  COMPUTER  AND  ELECTRICAL  ENGINEERING 

* 

*  20  APRIL  1987 

★ 

*  COMPUTES  THE  COMMUTATION  ANGLE  AND  LOAD  RESPONSE 

*  FOR  A  SINGLE  PHASE  FULL-WAVE  RECTIFIER.  USER 

*  INPUTS  ARE  LOAD  AND  SOURCE  CHARACTERISTICS,  FIRING 

*  ANGLE  AND  TRANSFORMER  INDUCTANCE. 

* 

★A**************************************************** 


*  VARIABLE  DICTIONARY 

★ 

a***************************************************** 


* 


* 


* 

* 


★ 


* 


A  FIRING  ANGLE  IN  RADIANS 

ALPHA  FIRING  ANGLE  IN  DEGREES 

DELTA  TIME-ANGLE  INCREMENT  IN  CURRENT 

CALCULATION 

ESIGN  DETERMINES  ZERO  CROSSING 
IN  MU  ESTIMATION  ROUTINE 
EFACT  ACCURACY  FACTOR  IN  METHOD  OF 
BISECTION 

ENEW  VALUE  OF  FUNCTION  DRIVEN  TO  ZERO 
IN  MU  ESTIMATION 

EOLD  VALUE  OF  FUNCTION  DRIVEN  TO  ZERO 

IN  MU  ESTIMATION 

ERROR  NEARNESS  TO  ZERO  REQUIREMENT  FOR 

FUNCTION  OF  MU  IN  BISECTION 
F  SOURCE  FREQUENCY  IN  HERTZ 

FX  SUBSCRIPTED  F  (IE.  FI,  ETC.). 

HOLDS  VALUES  GENERATED  IN  KX 
ROUTINE 

ID  LOAD  CURRENT  DURING  COMMUTATION 

IA  SCR#  1  CURRENT  DURING  COMMUTATION' 

I AC  SCR# 1  AND  LOAD  CURRENT  DURING 

CONDUCTION 

IAP  ARRAY  HOLDING  PHASE  A  CURRENT  VALUES 

FOR  ZETA  PLOT  ROUTINE 

IDP  ARRAY  HOLDING  LOAD  CURRENT  DURING  COMM’.  TATI.  N 

FOR  ZETA  PLOT  ROUTINE 
J  COUNTER 

KX  SINGLE  SUBSCRIPTED  K.  FUNCTION 

CONTAINING  CONSTANTS  OF  INTEGRATION 
IN  CURRENT  EQUATIONS 

KXX  DOUBLE  SUBSCRIPTED  K.  INTERMEDIATE 


1  7 


idt*  vrir*  *jr?jrjf  wj*  w  v  v  v  wv  w.vw  ^v.v.vr/^v^^.’VLV"/; 


-  *v  -  -  ¥-»  r.  -  i  <-*  w 


A 


★ 


★ 


i 


★ 


* 


RESULT 

L  LOAD  INDUCTANCE  IN  HENRIES 

LL  TRANSFORMER  INDUCTANCE  IN  HENRIES 

MU  COMMUTATION  ANGLE  IN  DEGREES 

MUST  START  ESTIMATE  OF  MU  FOR  BISECTION 
MUSTP  STOP  ESTIMATE  OF  MU  FOR  BISECTION 
MUOLD  INTERMEDIATE  RESULT  IN  BISECTION 

USED  IN  ACCURACY  CHECK 
MUNEW  INTERMEDIATE  RESULT  IN  BISECTION 

USED  IN  ACCURACY  CHECK 
PHI  LOAD  IMPEDANCE  ANGLE  IN  RADIANS 

PI  PI 

R  LOAD  RESISTANCE  IN  OHMS 

T1  TIME  CONSTANT  OF  COMMUTATION  CIRCUIT 

T2  TIME  CONSTANT  OF  CONDUCTION  CIRCUIT 

TEST  1  SIGN  DETERMINES  SIDE  OF  ZERO  CROSSING 
IN  BISECTION 

U  COMMUTATION  ANGLE  IN  RADIANS 

V  PEAK  VALUE  OF  SOURCE  IN  VOLTS 

VC  LOAD  EMF 

W  SOURCE  FREQUENCY  IN  R/S 

WT  INCREMENT  AI  ,LE  OF  HALF  CYCLE 

WTP  ARRAY  HOLDING  TIME-ANGLE  INFORMATION 

FOR  ZETA  PLOT  ROUTINE 
Z  LOAD  IMPEDANCE  MAGNITUDE  IN  OHMS 


» 


**  +  ***********■*★**★**********★*★*★★*****************★* 
* 


*  HOUSEKEEPING 

★ 


¥ 


****************************************************** 


1  REAL  L, LL, Kll, K12, K1 3 , Kl, K4  1 , K42, K4  3 , K4  4  , K4 

]  REAL  MU, K2, K3,MUST,MUSTP,MUOLD,MUNEW, ID, IA, IAC 

REAL  IA1, IA2, I A3 
REAL  IAP (182) , IDP (182) , WTP (182) 

1  ****************************************************** 
* 

*  INTERACTIVE  DATA  ENTRY 

* 

****************************************************** 


339  format;//) 

430  FORMAT (’ INPUT  DATA ( 1 )  OR  USE  INTERNAL ( 0 )’ ) 

format  (’data  may  be  input  in  real,  integer  i? 

(.EXPONENTIAL  (XX.  XE  XX)  FORM') 

■I  3  2  FORMAT  ('PEAK  SOURCE  VOLTAGE  (V)') 

‘•.9  3  FORMAT  ('  SOURCE  FREQUENCY  (HZ)') 

■'.  M  FORMAT  C  TRANSFORMER  INDUCTANCE  { H >  *  ) 

-l'.1)  FORMAT  (  1  LOAD  INDUCTANCE  ( H  )  '  ) 

4  3  6  FORMAT  ('LOAD  EMF,  VC  (V)') 

9  3  7  FORMAT ( ' LOAD  RESISTANCE  (OHM)') 

)  •<  FORMAT  ('  FIRING  ANGLE  (DEG)') 


18 


*  ^ > ’a  / 


'■■‘k  A .  JL 


\  A.  %  N 


i! 


FORMAT (' CONVERGENCE  FACTOR  (START  WITH  ABOUT  0.1)’) 
FORMAT ('RUN  IDENTIFICATION  NUMBER’) 

FORMAT ('200  ITERATIONS  AND  STILL  THE  FUNCTION  DOES 
;  NOT  CONVERGE ! ' ) 

FORMAT ('THE  ERROR  FACTOR  CURRENTLY  IN  USE  IS:  ’,G18.8) 

FORMAT ( ' THE  "ZERO"  YOU  ARE  TRYING  TO  MATCH  IS:  \G18.8) 

FORMAT ('THE  CLOSEST  YOU  HAVE  COME  SO  FAR  IS:  ’,G18.8) 

FORMAT (' INPUT  A  LARGER  ERROR  FACTOR') 

FORMAT ( 'MINIMUM  ALPHA  REQUIREMENT  VIOLATED.') 

FORMAT ( 'ALPHA  INCREASED  BY  1  DEGREE  TO :', F7 . 3 ,' DEGREES ' ) 


WRITE (9,400) 
ACCEPT  J 

I F ( J . EQ . 0 ) GOTO  10 
WRITE (9, 401) 

WRITE (9,399) 

WRITE (9, 402) 
ACCEPT  V 
WRITE (9, 399) 

WRITE (9, 403) 
ACCEPT  F 
WRITE (9,399) 

WRITE  (9,  404) 
ACCEPT  LL 
WRITE (9, 399) 

WRITE (9, 405) 
ACCEPT  L 
WRITE (9, 399) 

WRITE (9, 406) 
ACCEPT  VC 
WRITE (9, 399) 

WRITE  (9,407) 
ACCEPT  R 
WRITE (9, 399) 

WRITE (9, 408) 
ACCEPT  ALPHA 
WRITE (9,399) 

WRITE (9,409) 
ACCEPT  EFACT 
WRITE (9,399) 

WRITE ( 9, 398) 
ACCEPT  RUN 
GOTO  20 


NIT  I  A!  T7ATT~N  AN" 


•ENT  IF!  ('AT  I  CM  OF  "C  N..TAL 


*•****•***«*•••«*< 


F.FA  IT  -  0  .  1 


ALPHA=0 . 

F=60  . 

RUN=100 . 

20  WT=0 . 

1  =  1 

ENEW=0 . 

MUSTP=0 . 

J=0 

Tl=  (LL+2 . *L) / (2 . *R) 

T2=  (LL+L) /R 
PI=ACOS (-1 . ) 

W=2*PI*F 

Z=(R**2.+  ( (LL+L) **2.)  * ( W*  *2 . ) )  **.5 
PHI=ATAN2 (W* (LL+L) , R) 

30  A= (2.*PI/360.) * ALP HA 


ITERATIVE  CALCULATION  TO  ESTIMATE  THE  VALUE  OF  TH~ 
COMMUTATION  ANGLE 


10  F1=K1 (A,U,T1,W,V,LL, VC , R) 

F3=K3 (A,U,T1,W,V, VC, FI , PHI ,  Z ,  R) 
El=2 . *VC*T1 
E2=V/Z*SIN (A-PHI) 

E3=F3*EXP (- (PI-U) / (W*T2) ) 

E4=VC/R 

EOLD=ENEW 

ENEW=E1-F1-E2+E3-E4 
MUST=MUSTP 
MUSTP=U* 180 . /PI 
E=ENEW*EOLD 

ERROR=ABS (ENEW-EOLD) *EFACT 
IF(E.LT.O)  GOTO  100 
U=U+P 1/180 . 

GOTO  110 


METHOD  OF  BISECTION  TO  CALCULATE 
COMMUTATION  ANGLE 


V  V  V  '>  ">  ' 


vv,v:.v. 


V7\rif.'w  v/^ 


a 


WRITE (9, 440 ) 

WRITE (9, 399) 

WRITE (9, 441) EFACT 
WRITE (9, 399) 

WRITE (9, 442) ERROR 
WRITE (9, 399) 

WRITE (9, 443) ABS (E) 

WRITE (9, 399) 

WRITE  <9,  444) 

ACCEPT  EFACT 

ERROR=ABS (ENEW-EOLD) * EFACT 
ENDIF 
J=J  +  1 

U=MU*P I / 1 8  0 . 

F1=K1 ( A , U , T 1 , W , V , LL , VC , R ) 

F3=K3 (A, U, Tl, W, V, VC, FI, PHI, Z, R) 

E3  =  F3  *  EXP (- (PI-U) / (W*T2)  ) 

E=E1-F1-E2+E3-E4 

IF (ABS (E) . LE. ERROR) THEN 
GOTO  130 
ENDIF 

TEST1=E0LD*E 

IF  (TEST  1 . LT . 0 )  THEN 

MUNEW= (MU-MUOLD) / 2 . +MUOLD 
MUOLD=MU 
MU=MUNEW 
GOTO  120 
ENDIF 

MUNEW= (MU-MUOLD) *1 . 5+MUOLD 

MUOLD=MU 

MU=MUNEW 

GOTO  120 

1 30  U=MU*  P 1/18  0  . 

F 1  =  K 1  (A,  U,  Tl ,  W,  V,  LL,  VC,  R) 

F2=K2 ( A, T 1 , W, V, LL, VC, R, F 1 ) 

F3--K3  (A,  U,  Tl,  W,  V,  VC,  FI,  PHI,  Z,  R) 

,  *  *******************************************1 

CHECK  rCH  MINIMUM  FIR JUG  ATvGLE 

AC  H  EC  K  ~  2  *  V  *  3  I N  (A)  -  ( LL*  F 1 )  /Tl 

IF  (ACHECK  .  I.T  .  0  .  )  THEN 
ALPHA  =  AL1  HA+1 


•7»Oa 


I 


l>i' VXiMftjA  A  A  .>  lli  ,  ~v  ,■>  •,  . -- 

- 


WRITE (9, 445) 

WRITE (9, 446) ALPHA 

GOTO  30 

ENDIF 


CALL  OUTPUT (F 1 , F2 , E3 , T1 , T2 , V, VC, LL, L, R, ALPHA 
4,  MU, F, EF ACT , RUN) 

* 

»  CALCULATE  CURRENT  RESPONSE  THROUGH  COMMUTATION  AND 

*  ONE  CONDUCTION  PERIOD  ( WT=P I + ALPHA) 

* 

WT=A 

DELT  A=P 1/18  0 . 

140  IDP ( I) =F1 *EXP (- (WT-A) / (W*T1) ) -2 . *VC*T1 
I A I = V/ (W*LL) *COS (WT) 

IA2= (Fl/2)  *  EXP (- (WT-A) / ( W*T 1 )  ) 

IA3= (VC/LL) * (1 . -2 . *R*T1) * (WT/W) 

ZAP (I) =-IAl+IA2-IA3+F2 
WTP ( I ) =WT 
1  =  1  +  1 

WT=WT+DELTA 

IF (WT . LE . A+U) GOTO  140 

150  IAP ( I)  =  (V/Z)  *3  IN (WT-PHI ) +F3 *EXP ( - ( WT-A-U) / 

4  (W*T2 ) ) -VC/R 

WTP (I) =WT 
1  =  1  +  1 

WT=WT+DELTA 

IF  (WT.LE.  A  +  PDGOTO  150 

*  •  *  *  *  *  ****  ***-*********-**#***********-***#-**#**'******,***r 

* 

*  CREATE  PLOT  FILE  COMPATIBLE  WITH  ZETA  PLOTTER 

*  •  *  •  .  t******-****************************************** 

WTP  (16  1)  -A 
WTP ( 1  9  2  )  =  I  8  . 

N?T3=V /DELTA 
NPTS2-1 8  0  . 

CALL  PLOTS (0,0, L" PLOTS" ) 

CALL  PIGIN  ( 1 . 0, 4 . 0, 0) 

CALL  SCALE (IAP, 6. ,189,1) 

IDP (181) -IAP (191) 

IDP (182) -IAP ( 192 ) 

CALL  AXIS  (0 . , 0 . ,  'wt ' , -2, 10 . , 0 . , WTP (181), WTP (192)  ) 
'"ALL  AXIS  (0  0  'Cur  rent  (A)  '  ,  1 1 ,  6 . , 90 . ,  IAP  ( 1  8  1  ) 


->  u>  .v.v.M 


&,  IAP (182) ) 

CALL  LINE (WTP, IAP, NPTS2, 1, 0,2) 

CALL  LINE (WTP, IDP, NPTS , 1,0,2) 

CALL  PLOT (0 . , 0 . , 999) 

STOP 

END 

*  rt  ***  **ifcit->ti*‘*-*;fr*-*****:fr*-***'*****************1k*********** 

* 

*  DEFINED  FUNCTIONS  FOR  CONSTANTS  OF  INTEGRATION 

*  IN  MU  AND  CURRENT  EQUATIONS 

* 

REAL  FUNCTION  K1  ( A,  U,  Tl ,  W,  V,  LL,  VC,  R) 

REAL  Kll,  K12, K13, A,U, Tl, W, V, LL,  VC,  R 

K1 1=2 . / (1 ,+EXP (-U/ (W*T1) ) ) 

K12=V/ (W*LL) * (COS (A) -COS (A+U) ) 

K1 3=VC / (W*LL) *U* (1 . -2 . *R*T1) 

K1=K1 1  * (K12-K1 3+2 . *VC*T1 ) 

RETURN 

END 

REAL  FUNCTION  K2  (A, Tl , W, V, LL, VC, R, FI ) 

REAL  A, Tl , W, V, LL, VC, R,F1,K21,K22,K23 

K21=V/ (W*LL) *COS (A) 

K22  =  Fl/2  . 

K23= (VC* A) / (W*LL) *  (1 .-2 . *P*T1) 

K2=K21-K22+K23 

RETURN 

END 

REAL  FUNCTION  K3  ( A, U, Tl, W, V, VC, FI , PHI , Z, R) 

REAL  A,  U,  Tl,  W,  V,  VC,  FI,  K31,  K32,  K34  ,  Z,  R,  K33 

K3 1  =  F 1  *  EXP ( -U / (W*T1 )  ) 

K32=2 . *VC*T1 
K33=VC/R 

K3  4=V / Z  *  S IN (A+U-PHI) 

K3=K31-K32+K33-K34 

RETURN 

END 

*  •  *  *A********#****Hr-*-****i»-******'*#*****,it****i»*rt*****ifr** 

’  SUBROUTINE  TO  PRINT  INPUT  DATA  AND  CALCULATED 

■  CONSTANTS  TO  OUTPUT  FILE 

*  *  *  *  p  ***************************************************  * 


SUBROUTINE  OUTPUT (F 1 , F2 , F3 , T1 , T2 , V, VC, LL, L, R, ALPHA 
&  ,  MU,  F,  EFACT, RUN) 

.  0  FORMAT ( / ) 

>9  format;//) 

.  1  FORMAT (1  OX,  ' INPUT  DATA  &  CALCULATED  CONSTANTS 
&  PROGRAM  SINGLE  PHASE') 

2  FORMAT (3 4X, 'RUN  # ' , F 6 . 2 ) 

3  FORMAT (10X, 'PEAK  SOURCE  VOLTAGE  ',G18.8,' 

4  FORMAT (10X, 'SOURCE  FREQUENCY  ',G18.8,' 

5  FORMAT (10X, 'TRANSFORMER  INDUCTANCE  ',G18.8,' 

6  FORMAT (1  OX,  'FIRING  ANGLE  ',G18.8,' 

7  FORMAT (1  OX,  ' LOAD  INDUCTANCE  ',G18.8,' 

8  FORMAT ( 1  OX,  ' LOAD  RESTANCE  ',G18.8,' 

9  FORMAT (1  OX,  ' LOAD  EMF  ',G18.8,' 

0  FORMAT (10X, 'COMMUTATION  ANGLE  ',G18.8,' 

1  FORMAT (1  OX,  'ACCURACY  FACTOR  ’,G18.8) 

2  FORMAT (1  OX,  'COMMUTATION  TAU  ',G18.8,' 

3  FORMAT (1  OX,  'CONDUCTION  TAU  ',G18.8,' 

4  FORMAT (1  OX,  'CONSTANTS  IN  CURRENT  EQUATIONS') 

5  FORMAT (15X,  'Kl=' ,G18 . 8) 

! 6  FORMAT (15X,  ' K2=' ,G18 . 8) 


, G18 .8, 
, G18 .8, 
, G18 . 8 , 
, G18 .8, 
, G18 .8, 
, G1 8 .8, 
,G18 .8, 
, G18 .8, 
, G18 . 8 ) 
, G18 .8, 
, G18 .8, 


FORMAT ( 15X, 


FORMAT (15X, 'K3=' ,G18.8) 


WRITE  (7, 
WRITE  (7, 
WRITE (7, 
WRITE (7, 
WRITE  (7, 
WRITE (7, 
WRITE (7 , 
WRITE  (7, 
WRITE (7, 
WRITE (7, 
WRITE (7, 
WRITE (7, 
WRITE  (7, 
WRITE  (7, 
WRITE  (7, 
WRITE  (7, 
WRITE (7, 
WRITE  (7, 
WRITE (7, 
WRITE (7, 
WRITE (7, 
WRITE  (7, 
WRITE  (7, 
WRITE  C7, 
WRITE (7, 
WRITE  (7, 
WRITE  (7, 
WRITE (7, 
WRITE  (7, 
WRITE  (7 , 
WRITE  (7, 
WRITE (7, 


410) 

411) 

412)  RUN 
399) 

413)  V 
410) 

414)  F 
410) 

416)  ALPHA 
410) 

415)  LL 
410) 

417)  L 
410) 

418)  R 
410) 

4 1 9 )  VC 
410) 

422)  T1 
4  10) 

423)  T2 
410) 

420)  MU 
4  10) 


421) EFACT 
399) 

424) 

410) 

425)  FI 
410) 

426)  F2 
410) 


VOLTS  '  ) 
HERTZ ' ) 
HENRIES  '  ) 
DEGREES ' ) 
HENRIES ' ) 
OHMS  '  ) 
VOLTS ' ) 
DEGREES  '  ) 


SECONDS ' ) 
SECONDS  '  ) 
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INSTRUCTIONS  FOR  RUNNING  PROGRAM  SGLFAZ  ON  THE  UT  CDC  DUAL 

CYBER  170/750 


into  your  account,  then  type  the  commands  J  is  ted  he  low  at  the 
ropriate  prompt.  User  inputs  are  in  boldface  type. 

.READ  PF  ****  ****  SGLFAZ 

. FTN5  I=SGLFAZ 
. REWALLX 

. LDSET  LIB=ZETLIBF/ 

.  LGO 

,  ,,,,  ps  the  access  code  for  the  permanent  file  where  the  program  is 

V  ti  -J  _ 

this  point, you  will  be  presented  with  a  choice  of  using  input  data 
pained  in  the  source  file  or  entering  your  own.  If  you  choose  the 
ornal  data,  the  program  will  run  and  create  the  plot  and  output,  data 

y u  choose  to  input  your  own  data,  the  program  will  prompt  you  fro::: 
screen.  Correct  units  for  the  data  are  indicated  with  each  prompt. 

er  the  program  has  run,  the  load  response  can  bo  plotted.  At  the 
met  type: 

.DISPOSE  PLOT  ID=*  * 

is  the  unit  identifier  for  the  plotter  you  are  using. 

order  to  view  the  input  data  used  and  the  constants  calculate  •  in  the 
gram  (4,  kp,  k2,  kg,  a,  Ti ,  and  T2)  on  the  screen,  type  : 

.  SHOW  TAPE7 

you  want  to  print  the  data  on  a  line  printer  typo : 


.PRINT  TAPE7  ID  =  ** 
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